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Abstract. Shape Invariant potentials in the sense of [Gendenshte'in L.E., JETP Lett. 
38, (1983) 356] which depend on more than two parameters are not know to date. In 
[Cooper F., Ginocchio J.N. and Khare A., Phys. Rev. 36 D, (1987) 2458] was posed the 
problem of finding a class of Shape Invariant potentials which depend on n parameters 
transformed by translation, but it was not solved. We analyze the problem using some 
properties of the Riccati equation and we find the general solution. 



PACS numbers: 11.30.Pb, 03.65.Fd 
1. Introduction 

There has been much interest in the search of exactly solvable problems in Quantum 
Mechanics from the early days of the theory to date. To this respect, the Factorization 
Method introduced by Schrodinger [T^, [T7|, ^] and later developed by Infeld and Hull [|T^ 
has been shown to be very efficient. Later, the introduction of Supersymmetric Quantum 
Mechanics by Witten |f20| and the concept of Shape Invariance by Gendenshte'in [|ll| 



have renewed to great extent the interest in the subject. For an excellent review, see 

In particular, Shape Invariant problems have been shown to be exactly solvable, and 
it was observed that a number of known exactly solvable potentials belonged to such 
a class. The natural question which arose was whether all exactly solvable problems 
have the property of being Shape Invariant in the sense of |lTJ. This question has 



been treated in an interesting paper several years ago ||. There, the Natanzon class of 
potentials [HJ was investigated in detail. Following that line of reasoning, the authors 
gave a classification of Shape Invariant potentials whose parameters are transformed by 
translation. They proposed the general case which depends on an arbitrary but finite 
number n of parameters, and established the equations to be solved in order to find 
such a class. But they asserted to have failed to find any solution of the equations. 

f To appear in J. Phys. A: Math. Gen. (2000) 
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For several years this class of Shape Invariant potentials has been considered to 
be a good candidate to enlarge the class of known solutions of the Shape Invariance 
condition, see e.g. [Q, But the solutions are not known so far. 

As it seems to be an interesting problem we have analyzed it carefully and we have 
proved that it is possible to find in an easy way the solution. The main point is to use 
in an appropriate way some interesting properties of a related Riccati equation. As a 
consequence, the aim of this paper is to answer to the question proposed in |§. 

The organization of the paper is as follows. After a quick description of the problem 
of Shape Invariance in Section 2, we will develop in Section 3 the mathematical study 
of a particularly interesting first order ordinary differential equation system of key 
importance for the problem. Then we will proceed to study in Section 4 the problem 
of Shape Invariant potentials depending on n parameters. We will do some Ansdtze for 
the superpotentials assuming translations as the transformation law for the parameters, 
including the one proposed in || and its more immediate generalizations. The results 
are presented in some tables. 



2. Shape invariance and the Factorization Method 

We recall some basic ideas of the theory of related operators, the concept of partner 
potentials and Shape Invariance. Two Hamiltonians 

are said to be related whether there exists an operator A such that AH = HA, where 
A need not to be invertible. If we assume that 

A = ^- + W(x), (2) 

then, the relation AH = HA leads to 

W(V -V) = -W" -V , V -V = -2W' , (3) 
while the relation HA' = A^H leads to 

W(V -V) — W" — V' , V -V = -2W' . (4) 
One can easily integrate both pair of equations; we obtain 

V = W 2 -W + c, V = W 2 + W + d, 

where c and d are constants. But taking into account the equation V — V = —2W we 
have c = d. Therefore (see e.g two Hamiltonians H and H of the form ([!]) can be 
related by a first order differential operator A like (Q) if and only if there exists a real 
constant d such that W satisfies the pair of Riccati equations 

V — d — W 2 — W', V — d — W 2 + W' , (5) 

and then the Hamiltonians can be factorized as 

H = A*A + d, H = AA ] + d. (6) 
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Using equations in (|) we obtain the equivalent pair 

V-d= -(V-d) + 2W 2 , V = V + 2W'. (7) 

The potentials V and V are usually said to be partners. 

We would like to remark that these equations have an intimate relation with what 
it is currently known as Darboux transformations in the context of one dimensional or 
Supersymmetric Quantum Mechanics. In fact, it is easy to prove that the first of the 
equations @ can be transformed into a Schrodinger equation — 0" + (V(x) — d)<p = 
by means of the change — 0'/0 = W, and by means of 0'/0 = W the second of (|5|) 
transforms into —(f)" + (V(x) — d)0 = 0. The relation between V and V is given by fl7|). 
Obviously, 00 = 1, up to a non-vanishing constant factor. It is also worth noting that 
these Schrodinger equations express that and are respective eigenf unctions of the 
Hamiltonians (p]) for the eigenvalue d. These are the essential points of the mentioned 
Darboux transformations, as exposed e.g. in [13|, pp. 7,24]. 



The concept of Shape invariance introduced by Gendenshtein 11 : V is assumed 



to depend on certain set of parameters and equations @ define V and V in terms of a 
superpotential W. The condition for a partner V to be of the same form as V but for a 
different choice of the values of the parameters involved in V, is called Shape Invariance 
condition [|Il[ . 

More explicitly, if V — V(x,a) and V = V(x, a), where a denotes a set of 
parameters, Gendenshtein [|llj showed that if we assume the further relation between 
V(x, a) and V(x, a) given by 

V(x,a) = V(x,f(a))+R(f{a)), (8) 

where / is a transformation of the set of parameters a and R(f(a)) is a remainder not 
depending on x, then the complete spectra of the Hamiltonians H and H can be found 
easily. Just writing the a-dependence the equations (H) become 

V(x, a) - d = W 2 - W , V(x, a) — d = W 2 + W . (9) 

Therefore, we will assume that V(x, a) and V(x, a) are obtained from a 
superpotential function W(x, a) by means of 

V(x,a)-d = W 2 (x,a)-W'(x,a) , V(x, a)- d = W 2 (x, a)+W'{x, a) . (10) 

The Shape Invariance property in the sense of |TJ] requires the further condition (|8]) to 
be satisfied. 

The relationship of a slight generalization of the Factorization Method developed 
by Infeld and Hull [Tj|] with the Shape Invariance theory has been explicitly established 
in 0. There, the following identifications between the symbols used in the Factorization 
Method and those of Shape Invariance problems were found: 

V(x, a) — d = — r(x, f(a)) — L(a) , (11) 

V(x, a) -d= - r(x, a) - L(a) , (12) 

W(x,a) =k(x,a), (13) 

R(f(a)) =L(f(a))-L(a). (14) 
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3. General solution of equations y 2 + y' = a, zy + z' = b 

We will study next the general solution of a certain first order ordinary differential 
equation system. It will play a key role in the derivation of the main subject in this 
paper. The system is 

y 2 + y'= a , (15) 
yz + z' — b , (16) 

where a and b are real constants and the prime denotes derivative respect to x. The 
equation ( |15"|) is a Riccati equation with constant coefficients, meanwhile (|16D is an 
inhomogeneous linear first order differential equation for z, provided the function y is 
known. The general solution of (jlB]) is easily obtained once we know the solutions of 



(|15|), e.g. by means of 



_ bf x exp{fy(ri)dri}dZ + D 

{) ~ expiree} ' (17j 

where D is an integration constant J7|. 
The general Riccati equation 
dy 

— = a 2 (x)y 2 + ai{x)y + a (x) , (18) 

where a 2 {x), a\(x) and Oq(x) are differentiable functions of the independent variable 
x, has very interesting properties. It is to be remarked that in the most general case 
there is no way of writing the general solution by using some quadratures, but one 
can integrate it completely if one particular solution yi(x) of ( fT8|) is known. Then, the 



change of variable (see e.g. [10, 14 



1 1 

u = , with inverse y = yi , (19) 

yi-y u 



transforms (|i~8D into the inhomogeneous first order linear equation 



du 

— = -{2a 2 yi + a 1 )u + a 2 , (20) 
dx 

which can be integrated by two quadratures. An alternative change of variable was also 
proposed recently 0: 

yyi .,, . uyi 

u = , with inverse y = . (21) 

Vi-y u + yx 



+ a 1 )u + a , (22) 



This change transforms (fTS] ) into the inhomogeneous first order linear equation 
du (2a 
dx V y 1 

which is integrable by two quadratures, as well. We also remark that the general Riccati 
equation (|T8| ) admits the identically vanishing function as a solution if and only if 
ao(x) = for all x in the domain of the solution. 
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But the most important property of Riccati equation is that when three particular 
solutions of (|IBD , yi(x), y 2 (^), Dz{x) are known, the general solution y can be 
automatically written, by means of the formula 

2/2(2/3 — 2/1) + 2/1(2/2 — 2/3) , OQ , 

V = — 7 ri— , ( 23 ) 

(2/3 -yi)k + y 2 - 2/3 

where A; is a constant determining each solution. As an example, it is easy to check that 
y\k=o = 2/1 > y\k=i = 2/3 an d that the solution y 2 is obtained as the limit of k going to 00. 
For more information on geometric and group theoretic aspects of Riccati equation see 
e.g. ililH. 



We are interested here in the simpler case of the Riccati equation with constant 
coefficients (15l). The general equation of this type is 



dy 2 , . 

— = a 2 y + a x y + a , (24) 
dx 

where a 2 , at and ao are now real constants, a 2 7^ 0. This equation, unlike the general 
Riccati equation (]T8|) , is always integrable by quadratures, and the form of the solutions 
depends strongly on the sign of the discriminant A = a\ — Aa a 2 . This can be seen by 
separating the differential equation (23j) in the form 



a 2 y 2 



dy dy 

dx . 



aiy+ao « 2 ((?/+^) 2 -4) 



Integrating (24) in this way we obtain non-constant solutions. 

Looking for constant solutions of ([24]) amounts to solve an algebraic second degree 
equation. So, if A > there will be two different real constant solutions, when A = 
there is only one constant real solution and if A < we have no constant real solutions 
at all. 



These properties may be used for finding the general solution of (|i~5|). For this 



equation the discriminant A is just 4a. Then, if a > we can write a = c 2 , where c > 
is a real number. The non-constant particular solution 

yi(x) = ctanh(c(x — A)) , (25) 

where A is an arbitrary integration constant, is readily found by direct integration. In 
addition, there exist two different constant real solutions, 

y 2 (x) = c, y 3 (x) = -c. (26) 

The general solution obtained using formula (|23|) , is 



B sinh(c(x — A)) — cosh(c(x — A)) 
V ^ X) = ° B cosh(c(a; - A)) - sinh(c(x - A)) ' ^ ' 

where B = (2 — k)/k, k being the arbitrary constant in (|23|) . Substituting in ( |I7|) we 
obtain the general solution for z(x), 

= ~ C {B sinh(c(x - A)) - cosh(c(x - A))} + D 
B cosh(c(x — A)) — sinh(c(x — A)) 
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where D is a new integration constant. 

For the case a = 0, a particular solution is 

yi(x) = ^hi> (29) 



where A is an integration constant. If we apply the change of variable (piD with y 1 given 
by (|29|) , then (|T5|) with a = transforms into du/dx = 0. Then, the general solution 
for (|T5|) with a = is 

»<*> = l + fl£-A) - (30) 



with A and -B being arbitrary integration constants. Substituting in (|T7| ) we obtain the 
general solution for in this case, 

= 1 + B(x-A) ' (31) 

where D is a new integration constant. 

If now a = — c 2 < 0, where c > is a real number we find by direct integration the 
particular solution 

yi(x) = — ctan(c(x — A)) , (32) 

where A is an arbitrary integration constant. With either the change of variable or 
alternatively (|2T|), with yi(x) given by (|32"D we get the general solution of (|1^) for a > 

B sin(c(x - A)) + cos(c(a; - A)) 

y( x ) = - c n 71 JV\ ~T( Jv\ ' ( 33 ) 

B cos(c(a; — A)) — sm(c(x — A)) 

where B = cF, F arbitrary constant. Substituting in ([T7|) we obtain the general solution 
for z(x) in this case, 

-{B sin(c(x - A)) + cos(c(a: - A))} + g 
B cos(c(x — A)) — sm{c{x — A)) 
where D is a new integration constant. 

These solutions can be written in many mathematically equivalent ways. We have 
tried to give their simplest form and in such a way that the symmetry between the 
solutions for the case a > and a < were clearly recognized. Indeed, the general 
solution of (jlBD for a > can be transformed into that of the case a < by means of the 
formal changes c — >• ic, B — > iB and the identities sinh(ix) = zsin(x), cosh(ix) = cos(x). 
The results are summarized in Table |l[ 

Looking at the general solution of (]T5|) for a > 0, i.e. equation (^7]), one could be 
tempted to write it in the form of a logarithmic derivative, 

y(x) = — log \B cosh(c(x — A)) — sinh(c(x — A)) \ . 
dx 

This is equivalent except for B — > oo. In fact, if we want to calculate 
lim — log \B cosh(c(x — A)) — sinh(c(x — A)) \ 

B^oo dx 
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we cannot interchange the limit with the derivative, otherwise we would get a wrong 
result. But this limit for B is particularly important since when taking it in (|27|), we 
recover the particular solution (p5[). A similar thing happens in the general solutions 
(^0|) and fl33|). When taking the limit B — > oo we recover, respectively, the particular 
solutions ( pS| ) and fl32|), from which we have started. Both of fl30| ) and (|33| ) can be 
written in the form of a logarithmic derivative, but then the limit B — > oo could not be 
calculated properly. 



4. Shape Invariant potentials depending on an arbitrary number of 
parameters transformed by translation 

We will try now to generalize the class of possible factorizations considered in |L2], [7|. 
We analyze the possibility of introducing superpotentials depending on an arbitrary but 
finite number of parameters n which transforms by translation. This will give in turn 
the still unsolved problem proposed in ||. 

More explicitly, suppose that within the parameter space some of them transform 
according to 

/(aj) =a,i-ei, Vie T, (35) 

and the remainder according to 

.A",) «, • VjeT, (36) 

where ruf = {1, n}, and e, ^ for all i. Using a reparametrization, one can 
normalize each parameter in units of that is, we can introduce the new parameters 

m i = - ) Vz G T , and rrij = , V j G V , (37) 

for which the transformation law reads, with a slight abuse of the notation /, 

f(m i )=m i -l, Vi = l, ...,n. (38) 

Note that with these normalization, the initial values of each m; are defined by some 
value in the interval (0, 1] (mod Z). 

We will use the notation m — 1 for the n-tuple m— 1 = (mi — 1, ni2 — 1, . . . , m n — 1). 
The transformation law for the parameters ([38]) is just a particular case of a more general 
transformation considered in R]. As a corollary of a result proved there we have following 
one. The problem of finding the square integrable solutions of the equation 

dh 
dx 2 

according to the generalization of the Infeld and Hull Factorization Method treated in 
/0, Sec. 3], is equivalent to that of solving the discrete eigenvalue problem of Shape 
Invariant potentials in the sense of / fO|/ depending on the same n-tuple of parameters 



— + r(x,m)y + Xy = 0, (39) 



m = (mi, m2, . . . , m n ) which transform according to ftSBj). 
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In order to find solutions for these problems, we should find solutions of the 

difference-differential equation 

,2/ -\ ,2/ x dk(x,m+l) dk(x,m) T , . 

k 2 (x,m + l)-k 2 (x,m) + — — 7 - + — V — " = L(m)-L(m+l) , (40) 

dx dx 

where now m = (mi, m-i, ■ ■ ■ , m n ) denotes the set of parameters and m + l means 
771+1 = (mi + 1, 1712 + 1, • • • , m n + 1), and L(m) is some function to be determined, 
related to R(m) by R(m) = L(m) — L(m + 1). Equation ( |40"D is essentially equivalent 
to the Shape Invariance condition V(x, m) = V(x,m — 1) + R(m — 1) for problems 
defined by (|38|) 0. We would like to remark that (|40|) always has the trivial solution 
k(x,m) = h(m), for every arbitrary function h(m) of the parameters only. 

Our first assumption for the dependence of k(x, m) on x and m will be a 



generalization of the one used for the case of one parameter introduced in [12 



k(x, m) — ko(x) +mki(x) , (41) 
where k and ki are functions of x only. The generalization to n parameters is 

n 

k(x, m) = g (x) + ^ m i9i(%) ■ (42) 

i=l 

This form for k(x,m) is exactly the same as the one proposed in [Q, Eqs. (6.24)] taking 
into account ( |57| ) and (p8|), up to a slightly different notation. Substituting into (^) we 
obtain 

L(m) — L(m + 1) 

n n n n n 

= 2 5Z m ^+*^) +^H*^) + 2 (#o + 2o$>,:) . (43) 

j'=l i=l j'=l i=l i=l 

Since the coefficients of the powers of each have to be constant, we obtain the 
following first order differential equation system to be satisfied, 

n 

9j + 9jY, 9i = c ^ Vj G {1, n}, (44) 

i=l 
n 

9o + 9o 9i = c o , (45) 

i=l 

where q, i G {0, 1, . . . , n} are real constants. 

The solution of the system can be found by using barycentric coordinates for the <7j's, 
that is, the functions which separate the unknowns g^s in their mass-center coordinates 
and relative ones. Hence, we will make the following change of variables and use the 
notations 



9cm [X ) 

n . 
i=i 



-X>(s), (46) 
i=i 

1 n 

Vj(x) = g d (x) - g cm (x) = -(ng,j(x) - ^g^xjj , (47) 

i=i 

n 

n < * 



n . 

i=i 
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where j G {1, . . . , n}. Note that not all of the functions Vj are now linearly independent, 
but only n — 1 since YTj=i v j = 0- 

Taking the sum of equations ( fH|) we obtain that ng cm satisfies the Riccati equation 
with constant coefficients 

n 9cm + ( n 9cm) 2 = nc cm . 
On the other hand, we will consider the independent functions Vj(x), j G {2, . . . , n} to 
complete the system. Using equations ( f47D and (|44|) we find 



1 - 
v 'j = -( n 9j-22& 

i=l 

= -(g'j - g'i + g'j - 9*2 + ■ ■ ■ + g'j - g'j + ■ ■ ■ + 9j - g' n ) 

T^9 cm Cj C C m ; 

and we will take the corresponding equations from 2 to n. The system of equations 
and P5| ) is written in the new coordinates as 

^cm + ( n 9cm) 2 = nc cm , (49) 
v'j + ^n5- cm = Cj - c cm , V j G {2, . . . , n} , (50) 

#0 + govern = C , (51) 

and therefore the motion of the center of mass is decoupled from the other coordinates. 
But we already know the general solutions of equation (|49|) , which is nothing but the 
equation (|15D studied in the preceding section with the identification of y and a with 
ng cm and nc cm , respectively. Therefore the possible solutions depend on the sign of nc cm , 
that is, on the sign of the sum Y^h=i Ci °^ a ^ ^he cons t an t s appearing in equations fl44[) . 
Moreover, all the remaining equations (^) and (|5l|) are linear differential equations of 
the form (|T6|) , identifying z as Vj or g , and the constant b as Cj — c cm or c , respectively. 
The general solution of these equations is readily found once ng cm is known, by means of 
the formula (|T7|) adapted to each case. As a result the general solutions for the variables 
ngcm, Vj and g are directly found by just looking at Table [I] and making the proper 
substitutions. The results are shown in Table 0. 

Once the solutions of equations (^), fl50| ) and (^1|) are known it is easy to find the 
expressions for gi(x) and ^o(^) by reversing the change defined by (flop and (f47|). It is 
easy to prove that it is indeed invertible with inverse change given by 



SlW = gcm{x) 



£ «,(*), (52) 



=2 



9j(z) = gcm(x) + Vj(x) , Vj G {2, . . . , n} . (53) 

For each of the three families of solutions shown in Table |2], one can quickly find the 
corresponding functions gi(x), go(x), and hence the function k(x,m) according to fl42|). 
The results are shown in Table |3|. 

We can now calculate the corresponding Shape Invariant partner potentials by 
means of the formulas fllO|) , (|13|) and (0) adapted to this case. The results are shown 
in Table g 
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Let us comment on the solutions for the function k(x,m) in Table [5] and for the 
Shape Invariant potentials in Table f| we have just found. It is remarkable that the 
constants q, cq, of equations (£0|), (|45| ) appear always in the solutions by means of 
the combination Co + YH=i m i c i- On the other hand, D does not change under the 
transformation rrii — ► rrii — 1 since it depends only on differences of the m^'s. As 
D , D 2 , . . . , D n are arbitrary constants, D = D + Y^i=2 ^i( m i ~~ m i) can be regarded 
as an arbitrary constant as well. It is very easy to check that the functions k(x, m) satisfy 
indeed flloD , just taking into account that nc cm = Yli=i °i an d that when nc cm = C 2 , 
^™ =1 q/C = C, meanwhile ^2^ =1 Ci/C = —C when nc cm = —C 2 . Obviously, for the 
case nc cm = we have Y^=i c « = 0. As we have mentioned already, fllOD is essentially 
equivalent to the Shape Invariance condition V(x, m) = V(x,m — 1) + R{m — 1), but 
this last can be checked directly. In order to do it, it may be useful to recall several 
relations that the functions defined in Table ^| satisfy. When nc cm = C 2 we have 

f + = C(l - f + ) = C(B 2 - l)h 2 + , h' + = -Cf+h + , 
when nc cm = 0, 

f' = -Bf 2 , h' = -Bf ho + l, 
and finally when nc cm = —C 2 , 

f_ = C(l + f 2 ) = C(B 2 + l)h 2 _ , h'_ = Cf_h_ , 

where the prime means derivative respect to x. The arguments of the functions are the 
same as in the mentioned table and have been dropped out for simplicity. 

When we have only one parameter, that is, n = 1, one recovers the solutions for 
k(x, m) = ko(x)+mki(x) shown in the first column of J7], Table 6], and the corresponding 
Shape Invariant partner potentials of Table 7 in the same reference. 

For all cases in Table £|, the formal expression of R(m) is exactly the same, but 
either ^27=1 c * = nc cm have different sign or vanish. Let us consider now the problem 
of how to determine L{m) from R(m). The method does not provide the expression of 
L(m) but of L(m) — L(m + 1). In fact, there is a freedom in determining this function 
L(m). Fortunately, for the purposes of Quantum Mechanics the relevant function is 



R(m), from which the energy spectrum of Shape Invariant potentials in the sense of [I I 
is calculated 0. 

However, let us show how this underdetermination appears. Since 

n n 

R{m) = L(m) - L{m + 1) = 2(c + m iCi ) + a (54) 

i=i i=i 

is a polynomial in the n parameters rrii, and we have considered only polynomial 
functions of these quantities so far, L{m) should be also a polynomial. It is of degree 
two, otherwise a simple calculation would show that the coefficients of terms of degree 
3 or higher must vanish. So, we propose L(m) = Ylij=i r ij m i m j + Ym=i s i m i + ^ where 
Tij is symmetric, = r^. Therefore, there are |n(n + 1) + n + 1 constants to be 
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determined. Then, making use of the symmetry of in its indices we obtain 



L(m) - L{m + 1) = -2/] r^-m* - ^ r iS -/]sj. 

i,j=l 'J = l * = 1 

Comparing with (|5|]) we find the following conditions to be satisfied 

n n 

~^2 r ij = Vz G {1, . . . , n} , and - = 2c . 

i=i i=i 

The first of these equations expresses the problem of finding symmetric matrices of order 
n whose rows (or columns) sum n given numbers. That is, to solve a linear system of n 
equations with \n{n + 1) unknowns. For n > 1 the solutions determine an affine space 
of dimension \n{n + 1) — n — ^n(n — 1). Moreover, for n > 1 the second condition 
determine always an affine space of dimension n — 1. The well known case of of n = 1 
|12| . [7j] gives unique solution to both conditions. However, the constant t remains always 
undetermined. 

We will try to find now other generalizations of Shape Invariant potentials which 
depend on n parameters transformed by means of a translation. We should try a 
generalization using inverse powers of the parameters rrii, we know already that for 
the case n = 1 there appear at least three new families of solutions (see Table 6 in 0). 
So, we will try a solution of the following type, provided ^ 0, for all i, 

- fi(x) - 
k(x, m) = ) — h g {x) + V m^x) . (55) 

i=l 1 i=l 

Here, fi(x), gi(x) and go(x) are functions of x to be determined. Substituting into (|40|) 
we obtain, after a little algebra, 

L(m) - L(m + 1) = - > — J - 2g > — — 

f-r^ rriilmi + l)mj{mj + 1) ' ra^m* + 1) 

_ 2 mj9jfi + 2 fli/i + 2m i + 1 <Mi + 

mjfm, + 1) ^ mi + 1 r-f mjfmj + 1) dx 

where the dots represents the right hand side of (f43|). The coefficients of each of the 
different dependences on the parameters have to be constant. The term 



fifjjl + m, + nij) 
rriiimi + l)mj(mj + 1) 

involves a symmetric expression under the interchange of the indices i and j. As a 
consequence we obtain that fifj = Const, for all i, j. Since i and j run independently 
the only possibility is that fi = Const, for all i £ {1, . . . , n}. We will assume that at 
least one of the fi is different from zero, otherwise we would be in the already studied 
case. Then, the term 

n „ 



-j* rriiimi + 1) 



Shape Invariant potentials depending on n parameters 



12 



gives us go = Const, and the term which contains the derivatives of the /j's vanishes. 
The sum of the terms 

n j, n j, 

2 y- 9jfi 2 V m i9jfi 

is only zero for n — 1. Then, for n > 1 the first of them provides us J^r=i ft = Const, 
and the second one, gi = Const, for all i E {1, . . . , n}. This is just a particular case of 
the trivial solution. For n = 1, however, we obtain more solutions; is the case already 
discussed in [12], 0. It should be noted that, in general, 

2 9jfi 2 rrijgjfi _^ 2 V ^ (1 - ^ ] 

i,j=l i,j=l «,i=l 

as one could be tempted to write if one does not take care. Using the last equation as 
being true will lead to incorrect results. As a conclusion we obtain that the trial solution 
k(x,m) corresponding to that of the case n = 1 admits no non-trivial generalization to 
solutions of the type (|5"5]) . 

It can be shown that if we propose further generalizations to greater degree inverse 
powers of the parameters m^, the only solution is also a trivial one. For example, if we 
try a solution of type 

k(x, m) = J2 — + + g (x) + J2 m i9i (x) , (56) 

i,3=\ J 1=1 i=l 

where hij(x) = hji(x), the only possibility we will obtain is that all involved functions 
of x have to be constant. 



Now we try to generalize ( f42|) to higher positive powers. That is, we will try now 
a solution of type 

n n 

k(x, m) = g (x) + ^ m i9i{%) + 22 m i m j e ij( x ) ■ ( 57 ) 
i=l *,i=i 

Substituting into (pHi ) we obtain, after several calculations, 

n n 

L(m) - L(m + 1) = 4 mimjm k e i3 e k i + 4 ^ m^m^e^ + g k ) 

i,j,k,l=l i,j,k,l=l 
n / n . 

+ 2 ^ m i m j ( ( Gkl + 9l ^ + 

i,j=l \k,l=l 
n / n 

+ 4:^ m i e ij I ^2 ( ekl + 9l "> + 90 
i,j=l \k,l=l 

n / n , n 

+ 2 ^ mi I gfj (e^ + gj) + — ^(e ifc + ft) 

i=l V i,fc=l fe=l 

+ (e<i + ft) ( ^ (eifc + ft) + 2g j + ^- I J2 + ft) + 2ft ) • (58) 
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As in previous cases, the coefficients of each different type of dependence on the 
parameters have to be constant. Let us analyze the term of higher degree i.e. the 



first term on the right hand side of (58). Since it contains a completely symmetric sum 
in the parameters rrii, the dependence on the functions should also be completely 
symmetric in the corresponding indices. For that reason, we rewrite it as 

n . n 

4 2^ m i m j m keije k i = - ^ miinjm k (eije k i + e jk e a + e fci e j7 ) , 

i,j,k,l=l i,j,k,l=l 

from where it is found the necessary condition 

n 

^2(eije k i + e jk eu + e ki eji) = d ijk , k 6 {1, . . . , n} , 

i=i 

where dij k are completely symmetric in their three indices constants. The number of 
independent equations of this type is just the number of independent components of a 
completely symmetric in its three indices tensor, each one running from 1 to n. This 
number is \n{n + l)(n + 2). The number of independent variables is |n(n + 1) from 
the symmetry on the two indices. Then, the number of unknowns minus the number of 
equations is 

-n(n + 1) - -n(n + l)(n + 2) = —-(n - l)n(n + 1) . 
2 o o 

For n = 1 the system has the simple solution en = Const. For n > 1 the system 

is not compatible and has no solutions apart from the trivial one = Const, for all 

i,j. In either of these cases, it is very easy to deduce from the other terms in ( j5q) 

that all of the remaining functions have to be constant as well, provided that not all of 

the constants vanish. For higher positive power dependence on the parameters m,'s 

a similar result holds. In fact, let us suppose that the higher order term in our trial 

solution is of degree q, i ^m^m^ ■ ■ ■mi q T ilt _ iq , where T ilt _ >iq is a completely 

symmetric tensor in its indices. Then, is easy to prove that the higher order term 

appearing after substitution in (|40D is a sum whose general term is of degree 2q — 1 in 

the 7?Vs. It is completely symmetric under interchange of their indices, and appears the 

product of Ti u ,„ t i by itself but with one index summed. One then has to symmetrize 

the expression of the T's in order to obtain the independent equations to be satisfied, 

which is equal to the number of independent components of a completely symmetric 

tensor in its 2q — 1 indices. This number is (n + 2(q — l))\/(2q — l)!(n — 1)!. The number 

of independent unknowns is (n + q — l)\/q\(n — 1)!. So, the number of unknowns minus 

the one of equations is 

(n + q-l)\ (n + 2(g-l))! 



g!(n-l)! (2q - l)!(n - 1)! ' 
This number vanishes always for n — 1, which means that the problem is determined and 



we obtain that T\ 5 jl = Const., in agreement with |12[ p. 28]. If n > 1, one can easily 



check that for q > 1 that number is negative and hence there cannot be other solution 
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apart from the trivial solution 7^...^ = Const, for all zi, . . . , i q E {1, . . . , n}. From 
the terms of lower degree one should conclude that the only possibility is a particular 
case of the trivial solution. 



5. Conclusions and outlook 



Let us comment on the relevance of the more important result of this paper, that is, 
the fact that we have been able to solve the differential equation system ( fHf ) and (|45|) . 
That problem was posed, but not solved, in a often cited paper by Cooper, Ginocchio 
and Khare in Physical Review D || pp. 2471-2]. They use a slightly different notation 
but it can be identified their formulas (6.24) with our (|42|) and our procedure by an 



appropriate redefinition of the parameters taking into account (137|) and (j35|). However, 
they failed to find any solution to these equations (for n > 2), and believed that such a 
solution could hardly exist. 

The conclusion is conceptually of great importance. That is, it has been made 
clear that an arbitrary but finite number of parameters subject to transformation is 
not a limitation to the existence of Shape Invariant partner potentials, and hence, to 
the existence of exactly solvable problems in Quantum Mechanics. This leaves the door 
open to the possibility to pose and maybe solve further generalizations. We have as well 
the possibility of englobing particular cases of known Shape Invariant partner potentials 
spread over the extensive literature in the subject (see e.g. |J and references therein) 
into one simple but powerful scheme of classification. In this sense, we think the solution 
we have found here is very important as it completes the excellent work started in ||. 

Other conceptual point of great importance is that we have gained much more 
generality in the solution to the problem by a particularly simple but powerful idea. 
That is, to consider the general solution of the Riccati equation with constant coefficients 
which gives all subsequent solutions, rather than particular ones. For doing this it has 
been of great use the important properties of the Riccati equation. 

As a byproduct of our present results and that of @ it is not difficult to see that 
for n = 1 most of the solutions contained in || Sec. VI], later reproduced e.g. in |J, 
are directly related to some results of the classic paper |12 |, since they are solutions of 
essentially the same equations. 
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Table 1. General solutions of the equations ( |l5|) and (|l6|). A, B and D are integration 
constants. The constant B selects the particular solution of ( |l5| ) in each case. 



Sign of a 




z(x) 


a = c 2 > 


n B sinh(c(x — A)) — cosh(c(x — A)) 


£{B sinh(c(x-A))-cosh(c(x-A))} + D 


C B cosh(c(a;-A))-sinh(c(a;-A)) 


B cosh(c(x — A))~ sinh(c(x — A)) 


a = 


B 


b{f(x-A) 2 +x-A) + D 


1+B(x-A) 


1+B(x-A) 


a = -c 2 < 


B sin(c(x — A)J + cos(c(x — A)) 


§{B sin(c(x-A)) + cos(c(x-A))} + D 


C B cos(c{x-A))-sirL{c{x-A)) 


B cos(c(x — A)) — sin(c(a: — A)) 
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Table 2. General solutions for the differential equation system (|49|), (p0| ) and (pl|). A, 
-B, £>o and Dj are arbitrary constants. The constant B selects the particular solution 
of ( [49l ) for each sign of nc cm . 



Sign of nc c 



ng cm (x) 



Vj(x) for j £ {2, . . . , n} and goM 



nc cm = C 2 > 



Cf+(x,A,B,C) 



cj c cm f + ( x A B> Q + Djh + {x,A,B,C) 



%f+(x, A, B, C) + D h + (x, A, B, C) 



nCcm = 



Bf (x,A,B) 



)h (x,A,B)+D j f (x,A,B) 



c h (x,A,B) + D f (x,A,B) 



-Cf- (x,A,B,C) 



^-^f-{x,A,B,C)+D j h-(x,A,B,C) 



(x, A, B, C) + D h- (x, A, B, C) 



where 

f It A D — g sinh(C(x-A))-co B h(C(x-A)) 

J+y-L, ^ ^ t B cos h(C(a:-A))-sinh(C(a;-A)) 



Mx,A,B) = JTS ± 



x-A) 



h + (x, A, B,C) — B cosh ( C (x~A))-sinh(C(x-A)) 
J(x-A) 2 +x~A 



ho(x, A, B) 



1+B(x-A) 



f (t A R r\ - B sin(C(,x-A)) + cos(C(x~A)) 
J-(- L '- n '- D ' 1 -' ) — B cos(C(x- A))-s\n(C(x~ A)) 



h-(x,A,B,C) 



B cos(C(x-A))-sin(C(x-A)) 
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Table 3. General solutions for k{x, m) of the form (f42|). A, B are arbitrary constants. 
D denotes the combination Dq + Y^i=2 Di( m i ~ m i)i where Dq, Di are the same as in 
Table ^. The constant B selects the particular solution of (49) for each sign of nc cm . 



Sign of nc cm 


k(x,m) = g (x) + T,7=i 




nc cm = C 2 > 


^(co + £™=l m lCi )f + (x, A, B, C) - 


h B,C) 


nc cm = 


(c + £™ =1 m t a)h (x, A, B) + (D 


¥B^ m *)f Q {x,A,B) 


nc C m = -C 2 < 


^(c + E?=i m lCl )f-{x, A, B, C) - 


\- Dh-(x,A,B,C) 



where /+ = f+(x,A,B,C), f = fo(x, A, B), /_ = /- (x, A, B, C) 

h + = h + (x, A, B,C), h = h (x,A, B), h_ = h_(x,A,B,C) are defined as in Table g 
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Table 4. Shape-Invariant partner potentials which depend on n parameters 
transformed by traslation, when k(x, m) is of the form (|4^) and m = (mi, . . . , m n ). 
The Shape Invariance condition V(x, m) = V{x, m — 1) + R(m — 1) is satisfied in each 
case. A, B and D are arbitrary constants. 

Sign of nccm V(x,m) — d, V(x,m) — d and R(m) when k(x,m) = go(x) + E™=i m i9i( x ) 



+ (D 2 - (B 2 - 1)( C0 + E?=i »w))ft 2 
(C ° +S - im ' C - )2 4 + #(2(co + E? =1 miCi) - E?=i c s )/ + h + 



+ (D 2 + (B 2 - l)(co + Er=i "w))/^ 
R(m) = L(m) - Lira + 1) = 2(c + £™ =1 m iCi ) + E?=i c i 



(co + E"=i "W) 2 ^ + (5 + B^p)(D + B(£%^ + l))/ 2 
+2(c + E?=i m iCi )(D + B( Sr = im - + i))/ fto - (co + E?=i rmc) 

(co + E"=i m lCl ) 2 ft 2 + (5 + + B(£%^ - l))/ 2 

+2(c + E?=i m iCi )(D + B( E? = im - - §))Mo + (co + E"=i m * c 

B(m) = L(m) - L(m + 1) = 2(c + E™ =1 ro»ci) 



-C 2 < - (C ° + e|^' C ' )2 /- + §(2( c o + £?=! miCi) + E?=i Ci)/-h- 

+ (D 2 - (B 2 + 1)(cq + E?=i m lCl ))h 2 _ 

- (C0+ S:^ C ' )2 ^- + #( 2 ( C + ^?=! - E?=i <H)f-h. 

+ (D 2 + (B 2 + l)(co + E?=i ™«Ci))fc 2 

R(m) = L(m) - L(m + 1) = 2(c + E" = i m<Ci) + E?=i <=i 



where /+ = f+(x,A,B,C), f = f (x, A, B), /_ = /- (a;, A, B, C) 

h+ = h+(x, A, B,C), ho = ho(x, A, B), h- = h-(x,A,B,C) are defined as in Table 



